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ABSTRACT 



Q ■ In two previous papers a semi-analytical model was presented for the hierarchical 

clustering of halos via gravitational instability from peaks in a random Gaussian field of 
density fluctuations. This model is better founded than the extended Press-Schechter 
model, which is known to agree with numerical simulations and to make similar predic- 
tions. The specific merger rate, however, shows a significant departure at intermediate 
captured masses. The origin of this was suspected as being the rather crude approxi- 
mation used for the density of nested peaks. Here, we seek to verify this suspicion by 
implementing a more accurate expression for the latter quantity which accounts for the 
correlation among peaks. We confirm that the inclusion of the peak-peak correlation 
improves the specific merger rate, while the good behavior of the remaining quantities 
£p-( is preserved. 

o 

Subject headings: cosmology: theory - galaxies: clustering - galaxies: formation 



H ! 1. INTRODUCTION 

Today it is widely accepted that the formation of cosmic relaxed objects proceeds via grav- 
itational instability from small to large scales. This clustering process can be studied in detail, 
following accurately the dynamics of collapse by means of N-body simulations, and statistically, by 
approximate semi-analytical models which enable a swifter and more efficient exploration of the 
parameter space. 

The most complete semi-analytical models developed to date are those by Lacey & Cole (1993, 
1994, hereafter LC93 and LC94), and by Manrique & Salvador-Sole (1995, 1996, Papers I and II). 
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All the quantities in both models are derived from the statistics of the filtered random Gaussian 
density field, <5(x) = [p(x) — (p)]/(p), of cold matter encountered at an arbitrary initial epoch, ti, 
after recombination when fluctuations are still linear. Similarly, in each model the dynamics of 
dissipationless collapse is approximated by the simple spherical model. These two models differ 
however in the following aspects. 

Lacey & Cole's (LC) model is based on the Press & Schechter (1974, PS) heuristic prescription 
for the mass functions of halos, duly extended to obtain the conditional mass function of halos at 
a given epoch subject to having a larger mass at a later time (Bower 1991; Bond et al. 1991). 
iV-body simulations show that this model gives highly satisfactory predictions (LC94). In spite of 
this remarkable success, there remains the caveat that the seeds of halos are arbitrary points and 
their associated collapsing clouds undetermined (possibly disconnected) regions around them (cf. 
LC93). Besides, the formation and destruction of halos are arbitrarily and ambiguously defined, 
leading to poorly motivated formation and destruction times (see Paper II). 

The model developed by Manrique & Salvador-Sole uses the new CUSP (ConflUent System 
of Peak trajectories) formalism, which enables one to follow the filtering evolution of peaks in the 
density contrast vs. filtering scale (<5 vs. R) diagram at ti, assumed to trace the mass evolution 
of halos. This model is better founded than the LC model insofar as high peaks are good seeds 
of massive halos (Bond & Myers 1996) and do follow spherical collapse (Bernardeau 1994). On 
the other hand, it provides a natural distinction between accretion and merger which makes the 
formation and destruction of dark halos well-defined and uniquely determined. 

As shown in Papers I and II both models provide very similar results for the mass function, 
the mass increase rate, and the merger rate at large captured masses. However, notable departures 
were observed in this latter rate at intermediate captured masses. This discrepancy was suspected 
as being caused by the rather crude approximation used in the CUSP model for the density of 
nested peaks, in connection with the cloud-in-cloud correction. In the present paper we seek to 
verify this point. We improve the nesting correction by taking into account the correlation among 
peaks and investigate the effects this has on the CUSP model. In § 2 we review the foundations of 
the CUSP model developed in Papers I and II. In § 3 we derive a more accurate approximation for 
the density of nested peaks and implement it in § 4. Our results are summarized in § 5. Comoving 
lengths are assumed throughout the paper. 



2. THE CUSP MODEL 

Based on the spherical collapse model, the peak ansatz states that there is a correspondence 
between non-nested peaks of fixed linear overdensity 5 C in the density field smoothed on scale R at 
any arbitrary initial epoch ti and quasi-steady halos of mass M at a time t >ti- The mass M of the 
spherical collapsing cloud, or equivalently, of the final halo associated with a peak is an increasing 
function of the filtering scale R, whereas S c is a decreasing function of t. The mass function of halos 



at t, N(M, t)dM, is therefore simply given by the density of peaks with 5 C on scales R to R + dR, 
Npk(R, 5 c )dR, properly corrected for nesting and transformed to the variables M and t through the 
relations R(M) and 5 c (t). 

The approximate expression for the nesting correction used in Paper I was (see next section 
for justification) 

N(R,5)=N pk (R,5)- n^p-N(R',6)N pk (R,S\R',6)dR'. (1) 

Jr (p) 

In equation (El), (p) is the mean density of the universe, M(R')/{p) is the volume, to Oth order in 
5, of the collapsing cloud associated with a peak with S on scale R' , and N p k(R, 5\R' , 5) dR is the 
conditional density of peaks with 5 on scales R to R + dR subject to being located at a point with 
5 on scale R' > R. This equation is a Volterra-type integral equation of the second kind for the 
scale function N(R, 5) of non-nested peaks, which can be readily solved in the standard way from 
the known functions N p k(R,8) and N p k(R,5\R',5) derived in Paper I (eqs. [11] and [17]). 

The functions R{M) and S c (t) required to obtain N(M, t) from N(R, 5) are fixed by simple 
consistency arguments (see Paper I) 

In the above equation D(t) is the linear growth factor at the considered collapse time t and 8 ce {t) 
is the density contrast for collapse at t linearly extrapolated to t, both functions being dependent 
on the particular cosmogony assumed. Constants q and <5 c o = <5 C e(*o) are therefore the only free 
parameters of the model. For each pair of values of these parameters, the filtering process translates 
into a different mass evolution of halos, with only one couple of values (approximately) recovering 
the evolution encountered in real gravitational clustering. The correct values of these two parame- 
ters can be obtained by adjusting the mass function predicted by the model to the mass function 
resulting from iV-body simulations at some given time, for example at to (this automatically leads 
to a similarly good fit at any other time; see §4). 

The remaining statistical quantities predicted by the model rely on the definitions of accretion 
and merger, and halo formation and destruction. Since in hierarchical clustering halos grow con- 
stantly by capturing other halos, the fact that they are quasi-steady systems implicitly presumes 
that some of the captures are tiny enough for the steady state of the capturing systems not to be 
essentially altered. Notable captures yielding a transient departure from steadiness establish then 
the frontier between the initial and final quasi-steady systems. Accretion is therefore defined as 
any apparently continuous and derivable mass increase along the temporal series of quasi-steady 
systems subtending at each step the mass of those which precede them. In contrast, a merger is 
any appreciable discontinuity in the mass increase along such a temporal series causing a significant 
departure from the steadiness of the growing system at that step. A halo is said to survive as long 
as it evolves by accretion, whereas when it merges it is said to be destroyed. Note that when a halo 



is captured by one that is more massive it merges. However, the capturing halo may just accrete 
it if the relative captured mass is small enough. Only those events in which all the initial halos 
merge and are destroyed give rise to the formation of new halos. As shown in Paper II, the assumed 
correspondence between halos and non-nested peaks allows one to naturally identify the filtering 
processes at U tracing all the preceding processes, as well as to derive their respective rates and 
characteristic times (see §4). 



3. A MORE ACCURATE EXPRESSION FOR NESTING 

In the original version of the CUSP model the density distribution around peaks was approx- 
imated by top-hat spheres. In addition, the conditional density of peaks with 8 on scales R to 
R + dR, having 8' on R' > R and being located within a peak on that larger scale, was approxi- 
mated by the conditional density of peaks with 8 on scales R to R + dR simply having 8' on scale 
R' . Under these circumstances, the density of peaks with 8 on scales R to R + dR located within 
non-nested peaks with 8' on scales R' to R' + dR' (R' > R) takes the simple form 



M(R') 



N^ S \R -> B!,6-> 8') dR dR' = /^ ' N(R', 8') dR' N pk (R, 5\R', 8') dR, (3) 



Equation (||) leads, for 8' = 8, to the scale function of non-nested peaks (eq. |J|) and, for 8' = 8 — d8, 
to the density of peaks which become nested in an infinitesimal decrement of 8 (eq. |9| below), 
respectively tracing the mass function of halos and the density of halos merging into more massive 
ones. Any inaccuracy in this expression, therefore, translates into these two quantities and any other 
that is related. To improve equation (^), we would need to drop the top-hat approximation and take 
into account the peak-peak correlation. Approximate expressions for this latter quantity have been 
obtained by Bardeen et al. (1986; hereafter BBKS) and Regos & Szalay (1995). Unfortunately, 
these are not useful for our purposes; being only valid for large separations compared to the filtering 
scale, while what is needed here is the peak-peak correlation for separations up to R'. 

An alternative, more accurate, expression for the density of nested peaks can be obtained, 
however, from the conditional probability function P?(8,x,R,r) of finding the density contrast 8 
on scale R at a separation r from a peak with 8 and curvature x on that scale (averaged over 
the orientations of the second order derivative tensor) provided by BBKS. Using this conditional 
probability function we can write 

N^ S \R -^R',8^ 8') dR dR' = ^^ N(R', 8') dR' f dr 3r 2 f°° dx'P x ,(8\ R') 

\p) Jo Jo 

/oo 
dS P~ 5 {8', x', R', r)N pk (R, 8\R', 8) dR, (4) 

-oo 

with r in units of qR, and P X >(8',R') = N(R' ,x' ,8')/N(R',8') the curvature probability function 
for non- nested peaks with 8' on scales R' to R' + dR' . Expression (Q) takes into account the real 
density distribution around peaks and, since this includes the distribution of neighboring peaks 



(on every scale), it also accounts for the peak-peak correlation as desired. It is not exact however, 
because the conditional mean density of peaks with 6 on scales R to R + dR subject to having 5 
on scale R' at a distance r from a peak with 5' and x' on R' is approximated by N p i i .(R, 5\R' , 5) dR, 
the conditional mean density of peaks with 5 on scales R to R + dR subject to just having 5 on 
scale R'. 

In the Appendix we detail the calculations that allow us to write equation (H) in a compact 
form identical to equation (y), but for the new function N^ st (R,S\R' ,5') (eq. [ A13|| ) instead of 



the conditional density function Npk(R,5\R',5'). In Figure 1 we plot these two estimates of the 
density of nested peaks. As can be seen, they are very different for R' close to R while they 
asymptotically coincide for R' S> R- This is well understood. There is no correlation among 
peaks on very different scales. Thus the conditional mean density of peaks with 5 on scales R to 
R + dR subject to being located at a point with 5 on scale R' at a distance r from a peak with 5' 
and x 1 on the same scale is very well approximated, in this case, by the simple conditional mean 
density N pk (R,5\R' ,6) dR used in equation (|) leading to N^ st (R,5\R' ,5'). On the other hand, 
the spherical top-hat approximation is also very good in the limit for R' ^> R. Hence, not only 
is the new function N^ st (R,5\R' ,5') entirely accurate at this limit, but it also coincides with the 
former function Np^^R, 5\R' , 5'). 



4. EFFECTS ON THE CUSP MODEL 

All the quantities predicted by the new version of the CUSP model are readily obtained, 
following the same derivations as in Papers I and II, by simply replacing the function N p k(R, S\R', 5') 
by Np£ st (R,5\R',5'). This leads to the following expressions. 

The mass function of halos is 

N(M,t) = N(R,8 c )—, (5) 

with N(R, 5) the scale function of non-nested peaks at 5, solution of the Volterra equation 

N(R, 5) = N pk (R, 5) - jH ^p- N(R', 5) N^ st (R, S\R', 5)dR' . (6) 

In Figure 2 we compare this mass function with that obtained in the LC model. The parameters 
5 c q and q entering in the CUSP model were fixed by adjusting it to the PS mass function in the 
range 5 x 10 13 M Q < M < 10 15 M Q at z = 0, where the latter gives a very good fit to the results 
of iV-body simulations. Both functions are very similar for massive halos at any redshift. Only for 
small masses does the CUSP mass function predict slightly fewer halos than that of the PS. Note 
also that the large differences between the new and old approximations for nesting at R' ~ R (Fig. 
1) have no appreciable effect on the mass function: the original and new solutions being almost 
exactly superimposed. This indicates that the mass function is particularly insensitive, in the mass 
range of Figure 2, to the peak-peak correlation. 



The instantaneous mass accretion rate for halos of mass M at t is 



,(M,t) 



(x- 1 ) dM 
a 2 R ~dR 



dSc 
dt 



(7) 



with 02 the second order spectral moment (eq. |A6|| ), and (x _1 ) the average inverse curvature for 
the distribution function P X (S, R). There is no similar prediction in the LC model which does not 
differentiate accretion from merger. Only the total mass increase rates arising from the two models 
can be compared (see below). In Figure 3 we plot the accretion rates obtained from the original 
and new versions of the CUSP model. Once again, both solutions are almost fully superimposed 
showing that this rate is also insensitive to the peak-peak correlation. 

The instantaneous specific merger rate for halos of mass M at t per infinitesimal range of the 
resulting masses M' > M is 



l (M -» M',t) 



N m {R^R',5 c ) dR' 



N(R,5 C 



dM' 



d5 c 
dt 



(8) 



where N m (i? — ► R',S C ) dRdR'dS is the density of peaks with 5 C on scales R to R + dR becoming 
nested into non-nested peaks with 5f = 5 C — d5 on scales R' to R' + dR' 



N m (R^R',5) 



M(R') 



d 5f [N(R',5 f )N^ st (R,5\R',Sf)} 



0) 



To compare this specific merger rate with that yielded by the LC model we must take into account 
the different merger definitions adopted in the two models. In the LC model, merger is any mass 
capture experienced by halos. In contrast, in the CUSP model only notable captures are considered 
true mergers, while small captures contribute to accretion. As shown in Figure 4, the CUSP specific 
merger rate (using both the original and new approximations for the density of nested peaks) fully 
recovers that of the LC in the asymptotic regime at large AM/M = (M' — M)/M, but notably 
deviates for small AM/M. This departure is in the expected direction: for small mass captures 
the specific merger rate increases monotonically in the LC model, while it shows an abrupt cutoff 
in the CUSP model. 

From Figure 4 we also see that the original and new versions of the CUSP model now give 
markedly different solutions. They coincide at very large AM/M, owing to the similarity of the 
density of nested peaks at R' 3> R, but the old solution begins to decline, deviating from the LC 
solution, at AM/M much greater than unity. This implies that some smooth distinction between 
tiny and notable captures is already operating at AM/M 3> 1, which is meaningless. This is not the 
case for the new solution which keeps closer to that of the LC model until AM/M ~ 1, then rapidly 
falls off implying the existence of a sharp effective frontier between tiny and notable captures at 
AM/M <^ 1. Nonetheless, the new solution is not fully satisfactory: it becomes slightly negative 
for AM/M smaller than the effective threshold for merger, and the latter tends to be unreasonably 
small or large for extreme values of M. This shows the shortcomings of the new approximation 
used for the density of nested peaks. 



All the remaining quantities predicted by the CUSP model are determined by the three pre- 
ceding ones which consequently fix the behavior of the whole model. The instantaneous total mass 
increase rate for halos of mass M at t, r mass (M, t), is equal to the mass accretion rate (eq. J7|) plus 
the mass merger rate, defined as 

/■oo 

r™ ass (M,t)= AMr m (M^M',t)dM'. (10) 

Jm 

The instantaneous destruction rate (or global merger rate) of halos of mass M at t is 

roo 

r d (M,t)= r m (M -> M f ,t)dM', (11) 

Jm 



while their instantaneous formation rate take the form 

dS r 



f/ „,^ Nf(R,5 c 



dt 



(12) 



N(R,S C ) 
where N* (R, 5) dR d5 is the density of non-nested peaks appearing from 5 to 5 — d5, 

Nf(R,S)= d 5 N pk (R,5) + r ^p-N(R',5)d s/ N^ st (R,5'\R',5)\ 5 , =s dR' 

Jr \p) 

+ d R M pk (5, R) - J" ^p N(R', 6) d R M p n k est (S, R\5, R>) dR' (13) 



lR (P) 
+ ^lN(R,S) lim Ar$ st (6,R\6,R') 



A ^N(R,5) lim MT t! 
(p) R'-*R p 



(note the error in the expression quoted in Paper II), with J\f pk (5, R) and N pk est (5, R\S, R') defined 
for a fixed R instead of a fixed 5 (see Paper I). To see that expressions (O) and ( |T3| ) can be written 
in terms of the first three quantities above, one must simply substitute N(M,t), r% iass (M,t), and 



r 



(M,t) (eqs. ||, @, and |flT| ) into the conservation equation 

^f- = rf[M(t),t] - r d [M(t),t] - d M r a mass (M,t)\ M=M{t) , (14) 

for the number density of halos per unit mass along mean accretion tracks, M{t), solutions of the 
differential equation 

^-=r a mass [M(t),t]. (15) 

Finally, the distributions (and typical values) of formation and destruction times can be readily 
inferred from the expression giving the cumulative spatial number density of halos at U with masses 
from M(ti) to M{U) + SMfc) surviving until t f 

N sur (t f ) = N[M(ti),U]8M(ti) exp(- / ' r d [M(t'),t']dA, (16) 



with M(t) the mass at t of such halos calculated along their mean accretion tracks, and 



SM(U) = SM(tf) exp 



9m r mass(M, t )\M=M(t') & 



(17) 



the element of mass at t% evolving, through mean accretion tracks, into 5M(tf) at tt. Indeed, by 
taking t L = to and tj = t in equation ( |i~6| ) we obtain the distribution of destruction times of halos 
at to with masses between Mq and Mq + SMq, with 5Mq arbitrarily small, 

Mt) - -]y^) ^T = r ^ M{t ^ exp{-£r rf [M(t'),t'] dt'}, (18) 

while, by taking into account that the cumulative spatial number density of those halos at to which 
preexist at t < to, N pre (t), coincides with the cumulative spatial number density of halos at tj = t 
surviving until tf = to, equations (|l^) and (|l7]) lead to the distribution of formation times 

*,(t) = ^i^ ^p = rf[M(t),t] ex P {- jf* rf[M(t%t']dt>y (19) 



5. CONCLUSIONS 

We have investigated the effects that the inclusion of the spatial correlation among peaks has 
on the CUSP model. The predictions for the mass function of halos, their mass accretion rate, 
and their specific merger rate for large captured masses, are very similar to those arising from 
the original version of this model, which were already satisfactory. On the contrary, the specific 
merger rate at intermediate masses changes appreciably, now becoming more consistent with the 
corresponding prediction by the LC model. We, therefore, confirm that the poor behavior shown by 
this latter quantity in the original version of the CUSP model was not intrinsic to its foundations, 
but due to the rather crude approximation originally used for the density of nested peaks. A more 
accurate expression for the peak-peak correlation than that used here is, however, required for the 
predicted specific merger rate to be fully satisfactory. 

An interesting characteristic of the CUSP model is that it distinguishes between accretion and 
merger, which allows one to unambiguously define the formation and destruction of halos respec- 
tively as the last and next merger they experience. Mergers are those mass captures restructuring 
the systems, and hence, establishing their morphological properties until a new merger takes place. 
Such a distinction is therefore especially well suited for investigating the origin of the morphology 
of halos and, by extension, that of galaxies and galaxy clusters. Since the PS formalism involves 
much simpler expressions than those appearing in the peak theory, implementing that distinction 
in the former would have the added advantage of simplicity. The results of the present work suggest 
how it should be done. Such a modification of the LC model is used, in Salvador-Sole, Solanes, &; 
Manrique (1997), to study the origin of the correlation between halo concentration and halo mass 
recently found in cosmological iV-body simulations (Navarro, Frenk, & White 1996). 

The present work has been supported by the Direccion General de Investigacion Cientifica y 
Tecnica under contract PB96-0173. P.S. has benefited from a grant of the Swiss Science Foundation. 



APPENDIX: DENSITY OF NESTED PEAKS 



The same reasoning leading to equation (||) allows us to write its restriction to peaks with 
curvature x to x + dx on scale R, 5 on scales R to R + dR, and 5' on scales R' to R' + dR' 

M(R') 



N^ st (x, R^R',5^ 5') dx dR dR' 



(p) 



N(R',5')dR' / dr3r 2 / dx' P x ,{5' , R') 
Jo Jo 



dS P~ 5 (5', x', R', r) N pk (x, R, S\R f , 5) dx dR. (Al) 



In the preceding expression we have taken the volume of the collapsing clouds associated with peaks 
delimited by a top-hat window with radius proportional to the filtering scale R. 



Likewise, the density function of non-nested peaks with 5 on scales R to R + dR and curvature 
between x and x + dx satisfies the relation 

/•oo 

N(x,R,5) =N pk (x,R,5) 



ft 



N^ st (x, R^R',5^5) dR', 



(A2) 



similar to equation (y). From equations (Q), (Al), and (A2) we have 



P X (6,R) 



N pk {x,R,6) 
N(R,S) 



1 



B 



dR'M^lN(R',8) f dr3r 2 [°° dx' P x < (5, R') 

\p) Jo Jo 



dd Pg(6,x , R',r] 



N pk (x,R,6\R',d) 



(A3) 



SK "' ' N pk (x,R,6) 

Under the approximation P X (5,R) = N pk (x,R,5)/N pk (R,5), the curvature moments for this prob- 
ability function are accurate to within a few percent. Since the function N pk (x,R,8)/N pk (R,5) 
is bell-shaped, relatively symmetrical, and close to a Gaussian it can be further approximated by 
a normal distribution with identical mean and variance, x and (Ax) 2 . In taking this approxima- 
tion we will also extend the domain of integration over x' in equation (Al) to the whole range 
— oo < x' < oo. This does not introduce any appreciable error since x is considerably greater than 
3Ax for any reasonable values of 5 and R and any realistic spectrum. 

The distribution function Pg(5,x,R,r) is a normal distribution with mean and variance equal 
to (BBKS) 



5 x (r) 



7<5 (ip V V 



7 2 V7 



w- 



1-7 2 



7'0 + 



V 2 ip 



u^ 



and 



[A5 x (r)] 2 = a 2 Jl 



(1-7 2 ) 



i? + 2 7 V> + — / * 

u z 



3tfj' V 2 V> 



u 2 r 



7M 2 



(A4) 



(A5) 



In the above expressions, tp = £(r)/£(0), with r in units of qR and £(r) the mass correlation function 
on scale R, ip' = dip/dr, u 2 = (qR) 2 a2/cro, 7 = <T 2 /(a"oa"2), and 

/■oo U2(i+1) 

°HR)= -^-^P(k)eM-k 2 R 2 )dk, (A6) 

Jo 2tt z 



10 



with P(k) the power spectrum of the density fluctuations. 



Given the normal character of both P x > and Pg, the integration over x' on the right hand-side 
of equation (^|) can be performed analytically, the result being 



N^ st (R -> tf,6-> 5') dRdR' = ^ ' N(R', 5') dR' 

"1 POO 



PL POO 

/ dr3r 2 / d5P~ s (5',R',r)N pk (R,5\R',5)dR. (A7) 

JO J— oo 



The new function Pg(8, R, r), equal to the probability function of finding the value 5 at a separation 
r from a peak with 5 in the density field smoothed on scale R, is a normal distribution with mean 

1 2 _ + ^2~ ~ 1 2 7H- 

1 — 7 Z V 7 it / 1 — 7 V ir 

and variance 



^) = 7^(I + ^-SN + ^). (A8) 



[A5(r)] 2 = a 2 {l 



(1-7 2 



V> 2 + 2 7 V> + 



u 2 J u 2 



_ ( '3V' V 2 ^\ 2 3(V>') 
o 



i/ 2 / 7U 2 



(Ax) 2 / V 2 V 



+7T^N + ^) \. (A9) 



This enables us to perform the integration over 5 on the right hand-side of equation (A7) also 
analytically, obtaining 



JV3f st (R ^R',6^ 6') dR dR' = M . \ ' N(R', 6') dR dR' f 3r 2 N pk (R,S\R',6',r)dr, (A10) 

(p) Jo 

[u — e(r)v'(r)Y 



(P) 
where 



cxp 



N AR , S \g,g,r)dR= g| ;f;4f " l 2 '|-ffl J ^fldfl (All) 

is the conditional density of peaks with 5 on scales R to R + di? subject to being located at a 
separation r from a peak with 5' on scales R' to -R' + d-R'. In equation (All) R* = v3o"i/c>2, and 



iT[7(r),x*(r)] is defined just as i?[7,x*] in Paper I (eq.[17]), but for the r dependence in 7(7") and 
x*(r) introduced through 

e(r)= T-^n ' ^^^l 1 -^^} • (A12) 

Therefore, by defining 

N^(R, S\R', 5') = I' 3r 2 JV^U, <5|i?', 5', r) dr, (A13) 

equation ( |A10|) adopts the desired compact form similar to equation @ 

A^*(^ ->#,<&-► 5') dP d# = ^p- N(R', 5') dR' N^ S \R, MR', <*') dR- (A14) 
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Fig. 1. — Comparison between the old (thin solid lines) and new (thick solid lines) mean density of 
peaks with 5 on scale R nested within peaks with 5' = 5 = 5 c q on a larger scale R' = R + AR. The 
curves represent different values of R (corresponding to the same masses plotted in Fig. 4). The 
cosmogony assumed in all Figures to illustrate the behavior of the model is an Einstein-de Sitter 
universe with a n = —2 scale-free power spectrum of density fluctuations normalized to o"o(8/i _1 
Mpc) = .5. 



13 




12 
log (M/M ) 



Fig. 2. — Mass functions at two different epochs predicted by the CUSP model using the old (thin 
solid lines) and new (thick solid lines) approximations for the density of nested peaks and by the 
LC model (thick dashed lines). The curves corresponding to the old and new versions of the CUSP 
model are almost superimposed. 
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Fig. 3. — Mass accretion rates at z = predicted by the CUSP model using the old (thin dot-dashed 
line) and new (thick dot-dashed line) density of nested peaks; both curves almost fully superpose. 
Also plotted are the total mass increase rates predicted by the old (thin solid line) and new (thick 
solid line) version of the CUSP model compared with the prediction by the LC model (thick dashed 
line). 
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Fig. 4. — Instantaneous specific merger rates at z = of halos for M equal to 10 9 M (dotted 
lines), 10 11 M© (dashed lines), 10 13 M© (long dashed lines), and 10 15 M© (dot-long dashed line) as 
a function of the relative captured mass AM/M. (a) predictions by the CUSP model using the old 
(thin lines) and new (thick lines) approximations for the density of nested peaks, (b) these latter 
solutions (thick lines) are compared to the specific merger rates predicted by the LC model (thin 
lines) . 



